Correlations in emitters coupled to plasmonic waveguides 
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We report on quantum, classical, and total correlations in a set of distant quantum emitters coupled via their 
interaction with the plasmon modes of a one-dimensional waveguide driven by an external laser field. The 
coupling of the emitters with the plasmonic modes and its influence on the collective decay rate suggest that 
entanglement does not play a significant role in the qubit dynamics. Rather, discord is the quantity that matters 
and should be harnessed as a resource. 
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I. INTRODUCTION 

Entanglement is a cornerstone of quantum theory It 
has strong practical implications for futuristic quantum tech- 
nologies lIS HTOl and, consequently, it has been investigated in 
a wide variety of physical systems lfTTHT4l . For a thorough 
description of this phenomenon, a formal framework has been 
introduced over the past decade to quantify quantum lfT31l and 
classical lfT6l correlations in a physical system. The former, 
being of a more general character than entanglement, has also 
been hinted to play a role in quantum information process- 
ing through the so-called mixed state-based quantum comput- 
ing El EQ. 

Quantum emitters (e.g., single molecules 1 19, 20| or quan- 
tum dots B2H4231 ) coupled to surface plasmons of conduct- 
ing nanowires or waveguides have been seen as a promis- 
ing hardware for quantum information processing because of 
the strong coupling that can be achieved between the emit- 
ters and the plasmons ETl . The collective spontaneous de- 
cay and the plasmon-mediated emitter-emitter coupling have 
been investigated for a two-qubit (emitter) system close to a 
nanowire l24l l25l . as has the resonance energy transfer me- 
diated by different plasmonic nanowaveguides l26l . Entan- 
glement of two qubits, mediated by a V-groove plasmonic 
waveguide, has recently been theoretically studied 11271 l28l 
and quantified via the concurrence (29). We address this prob- 
lem within a more general framework: we calculate the entire 
spectrum for the entanglement, classical, quantum, and total 
correlations BOU321 for emitters coupled to plasmonic modes. 

We demonstrate that it is not entanglement but the quan- 
tum correlations (measured by the discord) that provide the 
relevant robust features that could be harnessed for plasmon- 
assisted information processing. Our results support the con- 
clusion that the entanglement of formation [29] is a natural 
metric and that, in the context of entropic measures, quantifi- 
cation via the concurrence might be inappropriate. 



II. QUANTUM EMITTERS, PLASMON MODES AND 
CORRELATIONS 

The Hamiltonian describing the dynamics of the two-qubit 
coupled interacting system [20] can be written in the compu- 



tational basis of direct product states \i) ® \ j) = 0, 1) as 
H$ = Hq + H\2, where Ho is the free particle term, and 

H 12 = \w(oP®oP + oP®gV), (1) 

gives the interaction energy, whose strength V depends of the 
configuration of the interacting dipoles. 

The dissipative dynamics is thus given by the total Hamil- 
tonian H = H$ + Hl, with Hi denoting the action of an exter- 
nal driving laser field, by means of the quantum master equa- 
tion (201 
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where o+ — |1,) (0,|, and a® = |0 ( ) are the raising and 
lowering Pauli operators acting on the jth emitter. Without 
loss of generality, we set Ta = T,- = F, and Ty = T*, = y(i^ 
j), the individual, and collective spontaneous emission rates, 
respectively. 



A. Collective decay effects 

The explicit form of the individual emitter decay rate T,, as 
well as the collective decay rate y associated to the dipole- 
dipole (qubit-qubit) interaction V depend on the particular 
physical setup under consideration. Although in what fol- 
lows we focus on the coupling of the emitters to the plasmonic 
modes of a channel waveguide, we start by briefly revising the 
main result for the interacting 'bare' quantum emitters. 



1. Bare quantum emitters 

We consider a pair of emitters coupled with interaction 
strength V, and individual transition frequencies in inter- 
action with a quantized radiation field in the vacuum state. 
The two emitters are separated by the vector r\2 and are char- 
acterized by transition dipole moments /I, = (0,-|D, 11,), with 



2 



dipole operators D,, and spontaneous emission rates T,. Here, 
the interaction strength and the incoherent decay rate read 
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where z = nfco'"i2> n denotes the medium refractive index, £q = 
coq/c, and coo = (coi +o>2)/2. 



2. Emitters coupled to a plasmonic waveguide 

We next consider that the emitters are located close to a 
plasmonic waveguide. The new qubit emission properties, as 
well as their dipolar interaction strength, depend on the cou- 
pling to the plasmonic modes. The dipole-dipole interaction 
strength and the incoherent decay rate can be calculated from 
the Green's function G(ri,r2) that describes the electromag- 
netic interaction between the two dipole moments (emitters) 
jUi 2 of frequency (0\ = <xh = (Oq as follows (see |25| for a 
discussion of the evaluation of these terms) 
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where the superscript "pi" indicates that the collective param- 
eters are now modified by the interaction with the plasmonic 
waveguide. As before, we set V^J = rf, and r P j = (r? 1 )* = 

y pl , with ij =1,2. 

It has been recently shown ll27l |28l that both the dipolar 
interaction and the collective damping for a pair of emitters 
close enough to a plasmonic waveguide hold simple analytical 
forms when the dominant contribution for emission is due to 
the plasmons. Thus, the Green function can be approximated 
by the plasmonic Green function contribution, G(ri,r2,&>) w 
G^ri.ia.co). 

If we set the plasmon wavelength Api = 2n/k p \, we 
have G5JE0 

V^ = ^^an(23cQ, f l = r pl j3 cos (2nQ , (5) 

where j3 = J3exp[— Api£/(2L)], with £ = d/Xpi, d is the dis- 
tance between the emitters, IT 1 = rf = Vf, L is the propa- 
gation length of the propagating mode, and the j3 factor that 
measures the fraction of emitted radiation by the propagating 
mode. It is worth mentioning that, due to radiative contri- 
butions, such a 'plasmonic approximation' breaks down for 
emitters separations shorter than ~ Api/4 [25]. 



We note from these two relations that | V pl |< r pl /2 and 
I 7P 1 I < TP 1 , which means that the interaction that arises due to 
the coupling to the plasmons is weak and the most important 
contribution due to collective effects comes from the damping 
yP 1 . 



B. Correlations 

The definitions involved in the calculation of the correla- 
tions are as follows. The quantum mutual information de- 
scribes the whole content of correlations in a given quantum 
system ||30H32l . It has been shown that quantum correlations 
(entanglement included) |[T5l and classical correlations lfl6l . 
in the sense of entropic measures, add up to give the quantum 
mutual information ||3T1 . Furthermore, this point has been re- 
cently emphasized, via the use of the relative entropy, within 
a unified framework that captures both quantum and classical 
correlations within the quantum mutual information ll32l . For 
a bipartite system, this can be written as: 



I(Pab) = S(pa)+S(Pb)-S(pab), 



(6) 



where S(p) = — Tr(plog 2 p) is the von Neumann entropy of 
density matrix p . 

A measure of classical correlations was introduced in 
Ref. Ifl6l as the maximum extractable classical information 
from a subsystem, say A, when a set of positive operator val- 
ued measures [7| has been performed on the other subsystem 
(say B): 



CC(p. 
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where S(p^) is the entropy associated to the density matrix 
of subsystem A after the measure. Such correlations must be 
non-increasing, and invariant under local unitary operations, 
and CC(p A b) — if and only if p AB ~Pa®Pb- 

If {|0) , |1)} define the basis states for the qubit B, the pro- 
jectors can be written as FFJ = 1 ® \j) j — a,b, where 

\a) = cos0|O) +e^sin0|l), \b) = sin0 10) -cos0|l), 
and the optimization is carried out over angles 9 and . The 
measure CC is antisymmetric by definition, and, without loss 
of generality, we take the qubit B to be the one measured. 

Following the definition for CC(pab), a simple way to de- 
fine the total quantum correlations in a bipartite system is 
D(Pab) = I(Pab) — CC(Pab). In terms of the von Neumann 
entropies, the quantum correlations, which coincide with the 
definition for the quantum discord given in Ref. [15], read 

D(pab) = S(p B ) - S(pab) + inf £>S(p A |n* ) • (8) 

{ n j} j 1 

For pure states, D = S(pb), and D = iff the system is purely 
classically correlated. 

Quantum entanglement is quantified by the entanglement 
of formation 



EoF(p) = h 



1 + V /1-C 2 (p)^ 



(9) 
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where h (x) — —x log, x — (1 — x) log 2 (1 — jc) denotes the binary 
entropy function ll29l . 

Additionally, consider, in decreasing order, the eigenvalues 
A,- of the matrix ^/PabPab, where p AB = (a y (gs cr v )p /1B (<7 y <g> 
(J v ) and Pab is the elementwise complex conjugate of p. The 
concurrence C can be defined as 



C(Pab) =max{0,Ai -A 2 -A 3 -A 4 }, 



(10) 



where the A/'s are as introduced above or, equivalently (also 
in decreasing order), the square root of the eigenvalues of the 
non-Hermitian matrix PabPab l29ll . 



III. DYNAMICS OF CORRELATIONS 

We calculate the exact dynamics of correlations by solving 
the master equation Q, and taking into account the defini- 
tions introduced in section MB I The notation for the correla- 
tions reported in the graphs throughout this paper is as fol- 
lows: total correlation or mutual information (purple, thick- 
solid line), classical correlation (green, doubly-dashed line), 
quantum discord (blue, thin-solid line), entanglement of for- 
mation (red, dashed line), and concurrence (pink, dotted line). 

To compare the physical scenarios described in Sec. [TTJ, we 
first present correlations when the emitters interact solely with 
the vacuum electromagnetic field, i.e., in the absence of cou- 
pling to plasmons. In terms of the collective effects, such a dy- 
namics depends on the separation between emitters and also 
on the orientation of their dipole moments. 

A typical situation where the dipole-dipole interaction is 
much greater than the incoherent part V ^> T ^ y [Fig.[TJa)] 
can be observed, for example, in a system of diluted molecules 
of Terrilene in a dispersive crystal |fl9l , for which a rich entan- 
glement dynamics has been predicted 11201 . Figure[TJa) shows 
that correlations reach high values and oscillate rapidly be- 
cause of the strong dipole-dipole interaction energy V; how- 
ever, the correlations decay within a short time and the system 
becomes uncorrected quite rapidly. 

A different scenario, depicted in Fig.[TJb), considers a pair 
of dipoles parallel to each other and perpendicular to their sep- 
aration vector, which allows the calculation of the collective 
parameters V and y directly from Eqs. ([3]). Although the cor- 
relations initially show smaller values than before, they persist 
for a longer time. 

An interesting issue that arises from Fig. [T] is the discrep- 
ancy between the entanglement of formation and the concur- 
rence as quantifiers of entanglement. Figure |TJb) shows that 
concurrence reaches much higher values (and decays much 
slower) than the EoF; more surprising is the fact that con- 
currence reaches, for almost any time, higher values than 
the total correlations, a result that contrasts with the defini- 
tion of the mutual information as a quantity that accounts for 
all the correlations (classical and quantum), entanglement in- 
cluded lf3T1 l32l . Thus, within this framework, concurrence 
can indicate results that are well above the EoF, and hence 
does not allow a direct comparison with other entropic mea- 
sures such as the quantum discord; in contrast, the latter can 
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FIG. 1. (Color online) Correlations for a two-emitter system in the 
absence of coupling to plasmons. (a) V = IT and y = 0.2r. (b) 
Dipoles parallel to each other and perpendicular to their separation 
vector; r\2 = f Ao, Ao = 2n/ko. The considered initial state is sep- 
arable: p(0) = 1 10) (10|. In all graphs presented in this work: Total 
correlation or mutual information (purple, thick-solid line), classi- 
cal correlation (green, double-dashed line), quantum discord (blue, 
thin-solid line), entanglement of formation (red, dashed line), and 
concurrence (pink, dotted line). 



be compared, on the same grounds, to the EoF 11331 . This is 
explicitly shown in Fig.[TJb). 

The influence of the plasmonic waveguide is illustrated in 
Fig. [2] and can be directly compared with the results in Fig.[T] 
According to the relations derived from Eqs. Q [25JE3, one 
of the parameters V pl or y? 1 can, in principle, be maximized 
by 'switching off the other one. For the interqubit distance 
d = | Api [Fig[2jb)], the plasmonic modes allow for an effec- 
tive enhancement of correlations, as can be seen from a di- 
rect comparison between the inset of Fig. |2|b) and Fig. |TJb). 
This said, note that the correlations decay more rapidly in the 
plasmon-assisted case because here the emitters distance is 
such that the collective rate is switched off, and the nonlocal 
effects are purely due to the weak dipole-dipole interaction 
V p . The correlations are enhanced as /3 tends to 1. For an 
emitters separation d = Api, the realistic value /3 = 0.82 can 
be obtained for L = 2/zm, and this is shown in Fig. |2ja): for 
this interqubit distance, V pl is switched off and the correla- 
tions reach larger values than before and are maintained for 
a much longer time. We also point out that for the parameter 
window considered in Fig.|2j classical correlations do not play 
a major role. 



FIG. 2. (Color online) Correlations for the emitters interacting with a 
plasmonic waveguide with the experimentally feasible set of param- 
eters j8 = 0.94, L = 2/lm, obtained for an operational wavelength 
A<) = 640 nm; this gives j8 ~ 0.82, and a plasmonic wavelength 
Api ~ 542 nm. (a) Separation between emitters d = Api, V pl = 0, 
(b) £ = d/Xp! = 3/4, yP 1 = 0, and I = 0.2 r pl . The inset corresponds 
to case (b), but in the absence of laser pumping. The notations for 
the graphs are as given in Fig.[TJ 



Interestingly, the quantum discord, which is in its essence 
different to entanglement, does play a role (and is more ro- 
bust than the EoF) during the qubit dissipative evolution. We 
remark that is not only the existence of quantum correlations, 
but the way they relate to each other, that matters, since the 
latter can lead to operational interpretations 11331 . In fact, 
Fig. [2^a) shows that for a time T = 10/r pl , entanglement, as 
calculated by the EoF, has almost vanished and the quantum 
discord quickly approximates to the value given by the total 
correlations. In contrast to this, Ref. G71 uses the concur- 
rence to show that a large amount of entanglement is gener- 
ated for this set of parameters. This can be seen in Fig. |2ja), 
where it is evident that the concurrence approximately equals 
the total correlations, in clear contrast to what is shown by the 
EoF. Within the formal framework of correlations presented 
here, the emitters entanglement is well below that reported in 
Ref. l27l . This is consistent with the approach of entropic 
measures followed here. 

Furthermore, it is straightforward to show that the den- 
sity matrix computed at t = T (Fig. 2|a)) corresponds to 
the mixed state p(r) = 0.9166 100) (00 + 0.0834 (¥" |, 
where = (|01) — 1 10) ) j \fl is the subradiant state. This 



FIG. 3. Optical control of correlations dynamics, (a) Total spectrum 
of correlations for laser amplitude i\ = —I2 = £ = 0.2 r pl , and initial 
state 1 1 ) . (b) Quantum discord dynamics for l\ = —£2 = £ = 0.4r pl , 
and initial state configurations ^/a|01) + \f\ — cc 1 10). (c) Compar- 
ison between QD, CC, and EoF for Of = 1. In all graphs, C = 1, 
V pl = 0, j8 = 0.94, and L = 2fim. 



state has a very small contribution from the entangled state 
l 1 ^^), and so it provides a small degree of entanglement. 
However, the mixture of this with the product state |0) <S> 
|0), produces a quantum correlated state with a type of 
correlation — the so-called quantum discord, which is differ- 
ent, in its very essence, to entanglement |[T5l[32ll . 

Although a larger value of j3 can enhance entanglement as 
well as the other correlations (not shown), we point out that 
care must be taken when choosing such a parameter; in par- 
ticular, the value j3 = 0.9 used in 11271 is not realistically at- 
tainable according to the plasmon dispersion relation reported 
there for d = Api. Indeed, if we include the 'best' achievable 
j3 (~ 0.94) factor, in the expression for j3 — > 0.9, the relation 
between the plasmonic wavelength and the propagation length 
should read Xpi/L ~ 0.08697, but from the known dispersion 
relation ll27l . no pair {Xpi,L} satisfies this constraint. This 
simple estimation leads us to conclude that, if j3 - — ' 0.94, the 
largest correct realistic value that should be used for the con- 
sidered emitters separation is ~ 0.82. 

The qubit correlation dynamics due to the collective effects 
arising from the emitters coupling to the plasmon modes, V pl 
and yf 1 , become even more interesting when they are pumped 
with a continuous laser field of amplitude £/ and frequency 
0), (targetting the emitter i). If we consider the laser excita- 
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tion to be in resonance with the emitters transition frequency 
(£>i = COq, stationary correlations can be obtained by making 
the amplitudes t\ = £2 = £, as shown in Fig. [2jb) for a dis- 
tance £ = 3/4, yP 1 = 0. In the same spirit, a higher stationary 
behaviour of correlations is obtained for £ = 1, by introducing 
a relative phase between the laser amplitudes: £\ = —£2 — £, 
as shown in Fig. [3ja). The correlations are much larger than 
those obtained in Fig.^b) because the laser excitation assists 
the very slow decay of the 'naturally created' antisymmetric 
state = (|01) - |10})/V2. This figure also shows that 
the concurrence reaches higher values than the total correla- 
tions, which might be seen as inappropriate when compared 
to the use of entropic metrics. 

We have shown that rather than entanglement, quantum dis- 
cord is the most robust correlation during the dissipative dy- 
namics, and, depending on the separation between emitters, 
the qubit-plasmon modes coupling can enhance the degree of 
correlations via the collective parameter "f l . For the sake of 
completeness, in Fig.[3jb) we consider, for a fixed laser ampli- 
tude, the role of the initial state preparation on the correlations 
dynamics. Here, we plot the quantum discord as a function 
of time, for initial states ^fa |01) + y/l - a |10), aeR. Al- 
though it is clear that, at t = 0, the discord takes the value 1 for 
a = 1/2, we have plotted it only for a range up to D = 0.3, in 
order to appreciate comparatively the correlations behaviour 
for all a. The trend in correlations is clearly influenced by 
a but, overall, the states converge to a common stationary 
value. The quantum discord, classical correlations, and EoF 
are plotted for the particular case a = 1 in Fig.[3jc), showing 
how an increment in the laser amplitude may produce a no- 
ticeable difference between the quantum discord and entan- 
glement (EoF). Thus, pumping with a continuous laser field 
can be used as an additional means to dynamically control the 
relationship between entanglement and the quantum correla- 
tions present in the emitters. Such a fact could be used as 
an important resource in the realisation of non-conventional 
quantum protocols ifTTl [T8l l34l . 

Bearing in mind that the quantum correlation, as measured 
by the discord, gives the amount of information that is not ac- 
cessible to local measurement, our results show that most of 
the information stored in the emitters is purely quantum (com- 
pare the thin-solid blue with the doubly-dashed green lines in 
all the figures). Also, we note that, in most cases, this type of 
information does not arise from entanglement, which is kept 
to a minimum, as can be seen in Fig.s[TJb),[2|a), and[3|c). This 
said, a possible optical control of the amount and class of in- 
duced quantum correlations can be carried out with external 
laser pumping, as depicted in Fig. [3] 



IV. DISCUSSION 



that the lifetime of the symmetric correlated state I 1 ?" 1- ) = 
(|01) + |10))/\/2 (an eigenstate of the bare Hamiltonian H$) 
is more robust (7 < 0). In the latter scenario, yP 1 = and 
the correlations are only mediated by the coherent interaction 
V p . The plasmon-mediated coupled emitters become more 
interesting because for large distances (d ~ Api), the incoher- 
ent decay yP 1 can take values close to T pl , holding the correla- 
tions for longer times. A particular case is shown in Fig.|2ja): 
d = Api, f 1 = 0.82IT 1 , hence VP 1 is 'switched off. The cor- 
relations are present even for times of the order of r pl f ^10 
thanks to the slow decay r pl — yP 1 of the antisymmetric state 
|<t>~). For this time scale, entanglement (as measured by the 
EoF) is almost zero, and the correlation that prevails is the 
quantum discord, which in turn tends to the same value as the 
total correlations. 

Aside from the physics that arises in the qubit-plasmon 
setup here considered, we highlight the discrepancy between 
the EoF and the concurrence as quantifiers of the emitters en- 
tanglement. This is due to the entropic origin of the EoF com- 
pared to that of concurrence (the difference between the two 
is explicitly stated in Eq.|9]l. Furthermore, in most of the plots 
of this work, the concurrence is larger than the mutual infor- 
mation (purple lines) in several different time frames. Since 
the quantum mutual information captures all the possible cor- 
relations (entanglement included), this result seems incorrect. 
Again, the explanation for this lies in the fact that the concur- 
rence does not have an entropic origin like the other correla- 
tion measurements calculated here. Strictly speaking, the con- 
currence is an entanglement monotone but is not an 'actual' 
entanglement metric in the sense that it obtains its meaning 
from its relation to the entanglement of formation and not the 
opposite l35l . With the emergence of quantum correlations 
beyond the entanglement, it is crucial to coherently quantify 
the latter in order to give a consistent physical interpretation of 
quantum phenomena. The understanding of whether a qubit 
system is entangled or discord-correlated is thus relevant for 
deciding the way in which the emitters could be operated as a 
physical device for performing information processing tasks. 

The relation between the correlations and the physical prop- 
erties of the hybrid system analyzed here is clearly showed in 
Figs. T||3 the lifetime of the correlations is enhanced thanks 
to the illumination with the coherent laser field; also, a large 
-factor (due to the plasmon channel) increases the degree of 
correlations, especially that of the quantum discord. It is ev- 
ident from Fig. [3] that an adequate manipulation of the light- 
matter interaction strength allows control of the dynamics fol- 
lowed by the quantum discord and entanglement, and there- 
fore of their metric value. This degree of quantum control is 
possible due to the interplay between the laser illumination 
intensity and the coherent emitter-emitter interaction l20l[36l . 



A comparison of the two physical scenarios considered here 
— emitters dissipative dynamics in the i) absence and ii) pres- 
ence of coupling to a plasmonic bus — leads us to the fol- 
lowing. In the former case, the correlations exhibit a longer 
lifetime in Fig. |TJb) because of the presence of the inco- 
herent interaction y for the distance r\2 = |Ao; this means 



V. CONCLUDING REMARKS 

We have computed entanglement, classical, quantum, and 
total correlations for a hybrid system composed of largely 
separated emitters coupled to the plasmonic modes of a one- 
dimensional waveguide, which is externally driven by a laser 
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field. We have illustrated, by direct calculation, that classi- 
cal correlations are the least assisted by the plasmon bus and 
that quantum discord is the dominant correlation that prevails 
throughout the whole dissipative dynamics and that this is en- 
hanced by the presence of the plasmonic collective excitations 
and by the external laser pumping. This tendency in the cor- 
relations has been analytically demonstrated for bare emitters 
coupled to the vacuum electromagnetic field in Ref . P6l . 

We have emphasized the entropic origin of the entangle- 
ment of formation, the quantum discord, and that of the quan- 
tum mutual information as quantifiers of correlations; we also 
computed the concurrence for direct comparison of our results 
with those reported in [27 28]. We found that the latter can 
reach values well above the EoF and the total correlations. 
This shows that care must be taken when confronting or using 
specific quantifiers of entanglement for describing a physical 
process, especially because the interpretation of the EoF as the 
cost of creation of an entangled state (with no regularization), 
leads to an upper bound on the actual degree of entanglement 
in a quantum system l37ll . 

We remark that a crucial feature of the mechanism devised 
here is the large separation between emitters (of the order of 
or greater than the plasmons wavelength), which could lead to 
long-distance quantum control and communication: the cor- 
relations are purely generated via the strong plasmon-emitters 
interaction, which in turn is reflected in the functional depen- 
dence of the collective parameters given by Eq. Q. An inter- 
esting by-product is the feasibility of tailoring the phase dif- 
ference that exists between the collective parameters, which 



might allow the switching on/off of the emitters correlations. 
Here, a parameter configuration where either V pl or goes to 
zero, an appropriate choice of initial conditions and laser tai- 
loring, enables the enhancement or suppression of the existing 
correlations. This result is beyond the scope of this work and 
will be addressed elsewhere ll38l . 

Recent experiments have demonstrated the control of light- 
matter interaction by means of plasmonic resonators 091 and 
silver nanowires j40l to enable the enhancement of the Pur- 
cell effect in quantum emitters, and the strength of the cou- 
pling to the plasmonic modes of such nanostructures. This 
offers the possibility of testing the results here reported with 
state-of-the-art technology. Furthermore, another experimen- 
tal demonstration of communication between two distant sin- 
gle emitters (organic molecules) via single photons has been 
recently reported lt4TTl . Thus, the setup proposed in this work 
has also the potential for demonstrating an additional degree 
of quantum control in which sensitive quantum information 
encoded in single photons can be transmitted and processed 
between coupled (but largely separated) single emitters. 
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